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Restricted open-shell Hartree-Fock (ROHF) theory is formulated as a projected self-consistent un- 
restricted HF (UHF) model by mathematically constraining spin density eigenvalues. The resulting 
constrained UHF (CUHF) wave function is identical to that obtained from Roothaan's effective Fock 
operator. Our a and /3 CUHF Fock operators are parameter-free and have canonical orbitals and or- 
bital energies that are physically meaningful as in UHF, except for eliminating spin contamination. 
The present approach removes ambiguities in ROHF orbital energies and the non- uniqueness of 
methods that build upon them. We present benchmarks to demonstrate CUHF physical correctness 
and good agreement with experimental results. 



Restricted open-shell Hartree-Fock (ROHF) theory 
was formulated by Roothaan some 50 years agoi A ma- 
jor drawback of this model is the lack of a unique effec- 
tive Fock operator^ Even though the ROHF wave func- 
tion and total energy obtained from different coupling 
schemes are the same, the resulting orbitals and orbital 
energies are different and lead to post-ROHF results that 
generally depend on them. The interpretation and phys- 
ical picture emerging from Roothaan's open-shell theory 
have always been somewhat blurry. Attempts to resolve 
these ambiguities, as well as many paradoxes resulting 
from them, are well documented in the literature.— ~— 

On the other hand, the physical picture of unrestricted 
HF (UHF) is clear J- It is a single-determinant wave func- 
tion with well-defined a and (3 orbital energies obeying 
Koopmans' theorem. It is straightforward to use it in 
post-UHF calculations by simply treating the a and (3 or- 
bitals explicitly and separately. The notorious problem in 
UHF, however, is spin contamination: the wave function 
is not an eigenfunction of S 2 . This weakness is ubiqui- 
tous and a serious detriment when bonds are stretched. If 
the UHF wave function suffers from severe spin contam- 
ination, as is the case when strong static correlation is 
present, then UHF is no longer a good starting reference 
point for post-UHF treatments of correlation or excited 
states. Once lost, good quantum numbers are hard to 
recover^ so when possible, it is preferrable to use ROHF 
as a starting point despite the ambiguities regarding its 
associated Fock operator. 

In recent work, we have developed a novel theory for 
treating strong correlations within an independent quasi- 
particle pictured— As a spin off of this work^ we have 
realized that the UHF energy can be written as a func- 
tional of the charge density matrix P = (~f a + 7^)/2 
and the spin density matrix M = (~f a — 7 /3 )/2, where 
■y a and 7^ are the a and (3 density matrices, respec- 
tively. Our proposed Constrained-Pairing Mean-Field 
Theory (CPMFT) departs from UHF by adopting a dif- 
ferent definition for M, a choice inspired by an underly- 
ing quasiparticle correlation picture. With this choice, 
CPMFT can accurately and efficiently describe static 
correlation and dissociate any molecule to fragments with 
ROHF energies, keeping the correct (S 2 ) all along the 



dissociation pathJ^ The connections between CPMFT, 
ROHF, and UHF turn out to be enlightening for formu- 
lating ROHF as a constrained UHF theory. The result- 
ing CUHF scheme here presented leads to well-defined a 
and (3 Fock operators with straightforward interpretation 
and no spin contamination. The ROHF wavefunction, 
energy, charge and spin densities remain the same; only 
the ROHF Fock operator is replaced by two UHF-like 
counterparts. As shown in benchmarks below, the mean- 
ing of the resulting orbitals and orbital energies is much 
more physical than in Roothaan's approach and provide 
a base for treatments of electron correlation and excited 
states. 

Theory. The familiar energy expression in ROHF is 
SROHP = 2^/^ i i+^/^ j (2a^iiHi)-^(ii|ji)), (1) 



where hij are one-electron integrals, (ij\kl) arc two- 
electron integrals in Dirac's notation, a and b are the 
coupling coefficients, and fi are the orbital occupations: 
1 for core (doubly-occupied, c) and for virtual (unoc- 
cupied, v) orbitals. In the case of high-spin open-shell 
systems under consideration, a = 1, b = 2, and / = 1/2 
for open-shells orbitals (singly-occupied, o). Roothaan's 
effective Fock operator is defined as 
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where F Q and are UHF a and (3 Fock matrices, and 
F cs _ (f« + f^)/2. At self-consistent field (SCF) conver- 
gence, all off-diagonal Frohf terms become zero. The 
choice of the diagonal elements in Eq.Q is completely 
arbitrary within a set of A and B coupling parameters: 
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Different values for these parameters have been suggested 
in the literature^ Although they do not affect the ROHF 
wave function and energy, they affect orbital energies 
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whose physical meaning is obscured because of this de- 
pendence. Choices guided to determine "canonical" sets 
that satisfy Koopmans' theorem may result in violations 
to the aufbau principle In this paper, these problems 
are resolved by abandoning the use of a single Fock op- 
erator. We will obtain the ROHF wave function by pro- 
jecting the UHF wave function self-consistently. Spin 
contamination in UHF is given byi£ 

S s = (S 2 ) - S Z (S Z + 1) = N - Tr( 7 V) (4) 

where S z = (N a — Np)/2 and N a (a = a, (3) is the num- 
ber of a electrons in the system. The previously proposed 
spin-constrained UHF (SUHF) approach^ introduces a 
Lagrange multiplier A in UHF to enforce S s = 0. How- 
ever, this is exact only in the limit of A — > oo. In this 
limit, the effective SUHF Fock matrices remain in the 
form of Eq. ([2]) We here propose an alternative method 
based on restricting natural occupations and spin density 
eigenvalues via finite Lagrange multipliers. 

In UHF, the natural occupations n are eigenvalues of 
P; they can be 0, 1, ^, or appear in "corresponding pairs" 
(n, I — This is a rigorous mathematical result follow- 
ing from P being the half sum of two idempotent density 
matrices^ In high-spin systems, the number of \ occu- 
pations is N a — Np = N s (we assume N a > Np always) . 
Note that TrP = (N a + Np)/2 = N e /2, where N e is the 
number of electrons. For clarity, we discuss below only 
the case where the number of orbitals TV is greater than 
N e but our results hold for N < N e too. The UHF j a 
are block-diagonal in the NO basis: 
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where N cp is the number of corresponding pairs and 
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and mi = \fn% — nf . The identity matrix in 7° accounts 
for unpaired electrons, traces to TV,, and is substituted by 
a corresponding zero matrix in ■y" . The other zero matrix 
has dimension N v = N — N s — 2N cp and corresponds to 
virtual (n = 0) unpaired orbitals. In the NO basis, M is 
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where, from Eq.©, M, = (7? - j?)/2 is 

I ° m n 

rrii 



(8) 



which is traceless with eigenvalues ±mj. The full spec- 
trum of M also includes | and eigenvalues, thus tracing 
to N s /2. Using the idempotency of j a and 7^, we get 



N, 



Tr(7 Q 7 ,3 ) = ^-2TrM^ 
Considering Eqs.([7]) and (J5J), it is evident that 



TrM 2 = 2]T; 



(9) 



(10) 



and hence 



Np - Tr(7 Q 7 /: 



(11) 



This readily means that to eliminate spin contamination 
in UHF all mi should be zero. Therefore, we propose 
to formulate ROHF as a constrained UHF scheme that 
enforces all m, to be zero. From Eq.©, m,; = implies 
that corresponding pair occupations become constrained 
to values of 1 and 0, thus effectively creating core (c) and 
virtual (v) orbital blocks. To enforce these constraints, 
we introduce Lagrange multipliers A^- and then write in 
a general basis 



-EblJHF — -E\jHF 



a Ma, 



(12) 



where the prime on the summation restricts it to cv and 
vc blocks. M is unconstrained in the 00 block and zero 
in other blocks. We next derive equations for Ay. 

The UHF energy is normally written as a functional of 
7" and 7' 9 . In our recent paper,— we have shown that 
the UHF energy expression can be alternatively written 
ms a functional of P and M, 

-Euhf = E cs + E c , (13a) 
E CS = 2J2 hijPij + ^(2(ij|fc0 - (ij\lk))P ik Pji (13b) 

ij ijkl 

E c = -^2(ij\lk)M ik M J i. 

ijkl 



(13c) 



E cs is the closed-shell energy expression given in terms 
of P, while E c is a "correlation energy" given in terms 
of M. The derivatives of £^ cs with respect to 7" and 7^ 
yield the usual closed-shell Fock matrix 
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On the other hand, the derivatives of E c are 
dE c , dE L _ 

kl 
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Hence, 
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which are the usual UHF Fock matrices. Now, the CUHF 
Fock matrices additionally require the derivatives of the 
constraints in Eq. (fl2|) with respect to 7" and 7^, which 
are trivially Ay/2 and —Ay/2, respectively. Defining 

A CUHF ag 



A CUHF 

ij 



Ag HF - ^ if {* e c a j e «}, 

or {iGdAjGc} 
otherwise 

yields the CUHF a and /3 Fock matrices, 
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The CUHF equations to solve are [F Q ,7"] = and 
[F^,7' 3 ] = 0. Subtracting these two SCF conditions and 
dividing it by 2 yields 

F C8 M - MF CS - A CUHF P + PA 



CUHF 



= 0. 



(19) 



Partitioning these matrices into core, open, and virtual 
blocks gives, 
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where we have used P cc = 1, P vv = P cv = P co = 
P t , D = 0, and P GO = M QO = |l. Together with Eq. (fT7|) . 
Eq. (|20cj) implies that A c „ = 2A C J J 1F at convergence. Dur- 
ing the iterative procedure, we choose this same value for 
A ct , because it guarantees 5 S = at each SCF cycle. Note 
that Eqs. ([20]) yield the SCF conditions for Roothaan's 
ROHF. Finally, our CUHF a and /3 Fock matrices are 
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These CUHF Fock matrices are different from the UHF 
ones only in the cv and vc blocks, and are different from 
Roothaan's effective Fock matrix of Eq.®. Our CUHF 
procedure yielding ROHF is surprisingly straightforward: 
one simply performs UHF with Fock matrices replaced 
by Eqs. ([2Tj) . These Fock matrices eliminate ambiguities 
arising in ROHF theory and produce a more physical 
UHF-like picture. In open-shell molecules, a and /? elec- 
trons feel different potentials; our F a and F 13 operators 
are different from each other and yield a orbitals dif- 
ferent from j3 orbitals that are true "canonical orbitals" 



TABLE I: Mean and mean absolute errors of ionization po- 
tentials (-ehomo in eV) of 24 open-shell systems. See the 
Supplementary Material for individual values.— 





ROHF" 


UHF 


CUHF 


ME 


7.38 


-0.68 


-0.54 


MAE 


7.38 


0.71 


0.61 



"Parameters from Ref. [21 



TABLE II: CN orbital energies in eV. In Roothaan's ROHF 
schemes (McWeeny and PGB), 5a is the open-shell orbital. 





3a 


4(7 


Itt 


5a 


6a 


McWeeny 


-33.81 


-16.89 


-13.87 


-6.21 


1.94 


PGB 


-33.46 


-16.44 


-13.68 


-18.17 


1.87 


CUHF a 


-34.54 


-19.87 


-14.26 


-16.06 


1.87 


CUHF P 


-33.46 


-16.44 


-13.68 


-1.88 


2.37 


UHF a 


-33.89 


-20.20 


-14.17 


-15.47 


1.87 


UHF /3 


-34.18 


-16.91 


-14.66 


-1.06 


2.42 


Exptl. IP a 




15.5 


14.4 


14.2 




"Ref.Q 



obtained by diagonalization. However, unlike UHF, they 
have no spin contamination, which is removed by La- 
grangian constraints. Their eigenvalues ef are physical 
orbital energies in the sense that they are associated with 
individual a and /3 orbitals, satisfy Koopmans' theorem, 
and the aufbau principle,— as opposed to many ROHF 
canonicalizations of Eq.([2]).— Our orbitals have previously 
been proposed in the literature as semi-canonical orbitals 
for MP2 and used in an ad hoc fashion^ Our present 
work shows that the Fock matrices for which these or- 
bitals are eigenfunctions appear from a constrained UHF 
optimization that eliminates spin contamination. 

Results. We have implemented CUHF in the Gaussian 
suite of programs^ and verified that our procedure 
converges to the ROHF energy. Unlike many ROHF 
schemes, CUHF has no issues with SCF convergence. 
This is undoubtedly related to the observance of the 
aufbau principle in our method. Since Koopmans' 
theorem is valid for CUHF, orbital energies approxi- 
mate ionization potentials (IP) and electron affinities 
(EA). In Table U we summarize the mean (ME) and 
mean absolute errors (MAE) of first IPs estimated via 
HOMO energies (ehomo) f° r 24 open-shell compounds 
selected from the G2 set.— Molecular geometries are 
optimized with B3LYP/6-31G(2df,p). CUHF results 
with a 6-311+-|-G(3df,3pd) basis are compared to UHF 
and the default (McWeeny) ROHF implementation^! in 
Gaussian. In all systems, the CUHF £homo captures the 
right physics yielding IPs comparable to those of UHF yet 
preserving the correct (S 2 ) expectation valued 

We have compared our CUHF orbital energies with 
those obtained by Eq.([2]) with parameters recently sug- 
gested by Plakhutin, Gorclik, and Brcslavskaya (PGB). 2 
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TABLE III: TDHF valence (V) and Rydberg (R) excitation 
energies (in eV) of open-shell molecules. Numbers in paren- 
theses are UHF spin contamination 5 S . 



System 


State 


UHF 


CUHF 


Exptl. 


BeF 


v 2 n 


4.20 


4.19 


4.14 


(0.001) 


R 2 E + 


6.34 


6.33 


6.16 




R 2 E + 


6.54 


6.54 


6.27 


BeH 


v 2 n 


2.69 


2.64 


2.48 


(0.002) 


r 2 n 


6.26 


6.25 


6.32 


CH3 


R 2 Ai 


6.54 


6.23 


5.73 


(0.012) 


R 2 A2 


7.73 


7.34 


7.44 


CO+ 


v 2 n 


6.93 


4.84 


3.26 


(0.141) 


V 2 E+ 


11.10 


9.81 


5.82 


CN 


v 2 n 


4.13 


0.95 


1.32 


(0.406) 


V 2 E+ 


5.42 


2.01 


3.22 


ME 
MAE 




1.43 
1.44 


0.45 
0.77 





The PGB parametrization is chosen to obey Koop- 
mans' theorem. However, the PGB scheme usually vi- 
olates the aufbau principle resulting in poor SCF con- 
vergence. Therefore, as a simple remedy, we have used 
the converged ROHF wave function and then diagonal- 
ized Eq.@ with PGB parameters in a single shot. The 
eigenvalues thus obtained are identical to those from the 
self-consistent PGB scheme. 

For the systems in Table U PGB gives the same first 
IP as CUHF, except for HCO whose ehomo are —10.88 



and —10.40 eV, respectively. The PGB scheme produces 
CUHF a virtual and /3 core orbital energies by construc- 
tion. However, the overall spectrum is appreciably differ- 
ent. In Table HH we present valence orbital energies for 
CN. In standard ROHF, 5er is predicted to be the open- 
shell orbital. Note the aufbau violation in PGB as previ- 
ously reported.— Results in Table iHl indicate that CUHF 
can well describe both a electron detachment and /? elec- 
tron attachment processes, yielding a spectrum that is 
fully interprctablc. 

Last, we present excitation energies of five small 
open-shell molecules calculated with time-dependent 
HF (TDHF) based on UHF and CUHF with a 6- 
311++G(3df,3pd) basis. The bond-lengths for BeF and 
CO+ (not included in the G2 set) are 1.355 and 1.078 
A, respectively. For TD-CUHF, we have used CUHF or- 
bitals and orbital energies in the TD-UHF procedure. 
Although this TD-CUHF scheme is not rigorous (one 
should perturb F, compute the response of P and include 
terms arising from A), this simple approximation turns 
out to be quite reasonable as shown in Table iHll When 
UHF spin contamination (5 S ) is small, TD-UHF and TD- 
CUHF give very similar results. As 8 S becomes larger, 
however, TD-UHF greatly overestimates the excitation 
energies. On the other hand, by retaining a spin pro- 
jected reference (S s = 0), TD-CUHF gives more reason- 
able excitation energies outperforming TD-UHF in spin 
contaminated situations. 

This work was supported by NSF(CHE-0807194) and 
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TABLE IV: £homo of open-shell systems (in eV). 



Atom 


ROHF 


UHF 


CUHF 


Exptl. IP 6 


H 


-3.40 


-13.60 


-13.60 


13.60 


Li 


-1.44 


-5.34 


-5.33 


5.39 


B 


-1.57 


-8.67 


-8.43 


8.30 


C 


-2.38 


-11.95 


-11.80 


11.26 


N 


-3.29 


-15.55 


-15.46 


14.54 





-4.87 


-14.21 


-14.37 


13.61 


F 


-6.55 


-18.54 


-18.62 


17.42 


Na 


-1.35 


-4.95 


-4.95 


5.14 


Al 


-1.22 


-5.94 


-5.72 


5.98 


Si 


-1.98 


-8.20 


-8.09 


8.15 


P 


-2.85 


-10.67 


-10.66 


10.49 


S 


-4.07 


-10.30 


-10.11 


10.36 


CI 


-5.40 


-13.09 


-13.00 


12.97 


OH 


-4.48 


-13.98 


-14.13 


13.01 


PH 2 


-2.89 


-10.25 


-9.94 


9.82 


SH 


-4.00 


-10.35 


-10.31 


10.37 


NH 


-3.25 


-13.82 


-13.79 


13.49 


O2 


-3.86 


-15.25 


-14.52 


12.07 




0.01 


1 n 4fi 




y.ou 


CH 3 


-2.01 


-10.46 


-10.18 


9.84 


C2H5 


-1.65 


-9.51 


-9.25 


8.12 


CN 


-6.21 


-14.17 


-13.68 


13.60 


HCO 


-2.60 


-10.73 


-10.40 


8.14 


CH30 


-3.93 


-12.16 


-12.29 


10.73 


ME 


7.38 


-0.68 


-0.54 




MAE 


7.38 


0.71 


0.61 





"Parameters from Ref. [2l| . 
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